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PERIODIC ORBITS IN HAMILTONIAN SYSTEMS WITH INVOLUTORY SYMMETRIES 


REEM ALOMAIR AND JAMES MONTALDI 


Abstract. We study the existence of families of periodic solutions in a neighbourhood of a sym¬ 
metric equilibrium point in two classes of Hamiltonian systems with involutory symmetries. In both 
classes, involutions reverse the sign of the Hamiltonian function. In the first class we study a Hamil¬ 
tonian system with a reversing involution R acting symplecticaUy. We first recover a result of Buzzi 
and Lamb showing that the equilibrium point is contained in a three dimensional conical subspace 
which consists of a two parameter family of periodic solutions with symmetry R and there may or 
may not exist two families of non-symmetric periodic solutions, depending on the coefficients of the 
Hamiltonian. In the second problem we study an equivariant Hamiltonian system with a symmetry 
S that acts anti-symplectically. Genetically, there is no S-symmetric solution in a neighbourhood 
of the equilibrium point. Moreover, we prove the existence of at least 2 and at most 12 families of 
non-symmetric periodic solutions. We conclude with a brief study of systems with both forms of 
symmetry, showing they have very similar structure to the system with symmetry R. 


1. Introduction 

A classical approach in the analysis of Hamiltonian systems is to study the existence of periodic 
orhits near equilibria. A basic theorem on the existence of periodic solutions in Hamiltonian sys¬ 
tems is the Liapunov centre theorem, which states that if the linearized flow at an equilibrium point 
has a simple purely imaginary eigenvalue satisfying a non-resonance condition then there exists 
a smooth 2-dimensional manifold which passes through the equilibrium point and consists of a 
one parameter family of periodic solutions, or nonlinear normal mode. In this work we extend this 
theorem to two classes of Hamiltonian systems with involutory symmetries where in both cases 
the involution reverses the sign of the Hamiltonian. In the first case, already studied by Buzzi and 
Lamb |5l, the involution is symplectic while in the second case it is anti-symplectic. 

In the literature, there are versions of the Liapunov centre theorem for reversible systems, but 
they mostly deal with the classical case where the reversing symmetry acts anti-symplectically. For 
example see Devaney (H . In this paper we consider two types of symmetry. First is the existence of 
periodic solutions in a time reversing Hamiltonian system equipped with an involution R that acts 
symplecticaUy. The problem was introduced and analysed by Buzzi and Lamb (5) . If the linear sys¬ 
tem has two pairs of purely imaginary eigenvalues, they prove in a neighbourhood of a symmetric 
equilibrium point the existence of a three dimensional subspace consists of a two parameter fam¬ 
ily of l?-periodic solutions with period close to 2n. In addition, they claim to find two families 
of non-symmetric periodic solutions whose period tends to 2n as they approach the equilibrium 
point, for an open dense set of coefficients (however there is a sign error in one of their calcula¬ 
tions). Motivated by this work, we looked at the problem using different coordinates, and hence a 
different set of invariants. We recover their result on the existence of symmetric periodic solutions 
but obtain a different conclusion for the non-symmetric solutions. We determine an expression in 
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the fourth order normal form and show that if this expression is positive there are two families of 
non-symmetric solutions, while if it is negative there are none. 

The second problem we discuss is the dynamics near an equilibrium point in an equivariant 
Hamiltonian system with an involutory (time preserving) symmetry S acting anti-symplectically. 
Bifurcations of equilibria in Hamiltonian systems with such symmetry have been considered re¬ 
cently by M. Bosschaert and H. Hanfimann (4| . Existence theorems for periodic solutions in sym¬ 
metric Hamiltonian systems can be found in Montaldi et al [14] , (TS) , but this and related work as¬ 
sumes the symmetry transformation acts symplectically. We prove that for systems with this anti- 
symplectic symmetry, generically, there are no symmetric periodic orbits in a neighbourhood of 
an equilibrium point. Moreover, we prove the existence of at least 2 and at most 12 non-symmetric 
families of periodic solutions (nonlinear normal modes) in a neighbourhood of the equilibrium 
point under the same generic conditions. 

In both cases, since the involution reverses the sign of the Hamiltonian and we assume the linear 
system is periodic, the equilibrium will be in 1:-1 resonance. 

The paper is organised as follows. In Section|2]we introduce basic facts and definitions of Hamil¬ 
tonian systems with symmetry. Section 3 lists normal forms of the Hamiltonian linear system L, the 
structure map / and the symmetry elements R and S in . Section 4 reviews the standard tool used 
to find periodic orbits in Hamiltonian systems: Liapunov-Schmidt reduction. In Section 5 we state 
and prove our theorem on the existence of families of periodic orbits in the f?-reversible Hamil¬ 
tonian system with R acting symplectically. In Section 6 we give our main result on the existence 
of periodic solutions in the S-equivariant Hamiltonian system with S acting anti-symplectically. 
Finally, in Section 7 we study the existence of periodic solutions in systems with the combined 
symmetry x Z| reversible/equivariant Hamiltonian system. 

2. Hamiltonian systems with symmetry 

In this section we recall some basic facts and definitions on Hamiltonian systems with sym¬ 
metry. Let iU^",a)) be a symplectic space, i.e an even dimensional vector space equipped with a 
symplectic form o). Recall that a symplectic form is a non-degenerate, skew symmetric, bilinear 
form. Then there exists a structure map / satisfying J* = -J [J* denotes the transpose of J) and 

= -I such that - {x,Jy) for x,y e where {.,.) is the standard inner product in 

Let H : — IR be a Hamiltonian function. The Hamiltonian vector field / generated by H is 

symplectic, i.e. its flow preserves the symplectic form o), and is defined by 


x = /(x)=/VH. (2.1) 

By using canonical coordinates for the symplectic form o) given in Darboux theorem (T] one can 
write 

0 -In 
\ In 0 , 

In this work we will deal with two types of symmetry, equivariant symmetries and time-reversing 
symmetries. 

Definition 2.1. Let S, R be two linear transformations of then 
(1) The vector field / is called S-equivariant if 

/(Sx)-S/(x),VxeR. 

If x(f) is a solution of II2.1L then Sx(f) is also a solution and S is referred to as a symmetry. 
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(2) The vector field / is called /^-reversible if 

fiRx)^-Rfix),yxeU. 

If x{t) is a solution of <2.1ll . then Rx[-1) is also a solution. Such a transformation is called a 
time reversing symmetry. 

The symmetry of a periodic solution is given by the following definition. 

Definition 2.2. Let x(f) be a periodic solution of the dynamical system x = fix). 

(1) If S is a symmetry of the system then x{t) is said to be S-symmetric if 

Sx[t + 6) = xit), 

for some 0 e S^. 

(2) If i? is a reversing symmetry of the system then xit) is said to be //-symmetric if 

Rxid- t) - xit), 

for some 0 e S^. 

Here we identify with IR/ TZ, where T is the period of xit). In both cases, a symmetric periodic 
orbit symmetric if and only if it is set-wise invariant. 

In the Hamiltonian context, (reversing) symmetries can arise in two ways: they can either be 
symplectic or antisymplectic. A (reversing) symmetry T is symplectic if(j)iTx,Ty) - (ji)ix,y),yx,y e 
R and anti-symplectic if a)iTx,Ty) = -coix, y), Vx, y e R. In matrix form we can choose a basis so 
that T is orthogonal, and then T is symplectic if TJ = JT and anti-symplectic if TJ = -JT. 

Note for example that by 12.IL if a reversing symmetry R is symplectic then it must reverse V//, 
and if we assume (as we may, and do) that //(O) = 0 then this is equivalent to HiRx) = -Hix), so 
that H is ‘anti-invariant’. There are in all 4 possibilities of symmetry, labelled as follows 


type 

(D 

/ 

H 

SE 

-1-1 

+ 1 

-1-1 

AR 

-1 

-1 

H-l 

SR 

-1-1 

-1 

-1 

AE 

-1 

+ 1 

-1 


Table 1. The ‘type’ refers to a transformation being Symplectic-Equivariant, or 
Antisymplectic-Reversing etc. 


Note that if T is an involution which reverses the sign of H, then any symmetric periodic orbit 
must lie in the set where H = 0. There may on the other hand be periodic orbits on which H is 
non-zero, and then T will exchange two such orbits, one with H>0 and the other with H <0. We 
will see this in more detail in later sections. 

3. Linear Hamiltonian systems with involutory symmetries 

In this section we give the normal forms of linear Hamiltonian systems with involutory sym¬ 
metries. Recall that an involution is a transformation of order 2. An important assumption that is 
required for studying the existence of periodic orbits is the presence of purely imaginary eigenval¬ 
ues of the linear Hamiltonian vector field. 

Let L e spjiZn, R) be a linear Hamiltonian vector field. Thus, 

LJ=-JL*, 
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where / is the structure map defined in the previous section. By Bochner’s theorem (3| , a (revers¬ 
ing) symmetry T can he chosen to he linear and orthogonal. Therefore, the (reversing) equivariant 
condition can he written as 

LT = ±TL, 

and the (anti-)symplectic property of T is given hy 

TJ=±JT. 


In (TT|, Hoveijn et al. gave normal forms of linear systems in eigenspaces of (anti-) automor¬ 
phisms of order two, which can he adapted to our problem. These normal forms are based on 
writing (/, T) - invariant subspaces. Since we are interested in generic systems with given symme¬ 
try, then by [TT] we can only focus on the case when L is semi-simple. Also, we assume that L has 
at least one pair of purely imaginary eigenvalues ±i. Normal forms of T,J and L are given in the 
following lemma. We use the notation. 


(1 0) 


0 -T 

( 0 1)’ 

J 2 - 

f 0, 


and S 2 - 


T 

0 



Lemma 3.1. LetL be a linear Hamiltonian vector field on 

i) SupposeL is R-reversible, withR acting symplectically (symmetry type SR). 

LetV be a minimal iL,J,R)-invariant subspace on which L has eigenvalues ±i. ThendimV = 
4 and RlY,Jly andL\y can take the following normal forms 


R\v - 






0 

h 



( J 2 0 ] 
0 J 2 j’ 


and L|y = 


/2 0 ) 

0 -/2 j ■ 


ii) Suppose now L is S equivariant, with S actinganti-symplectically (symmetry typeAE). 

LetV be a minimal iL,J,S)-invariant subspace on which L has eigenvalues ±i. ThendimV = 
4 and S\y,J\y andL\y can take the following normal forms 


f 0 S 2 'I 
S 2 0 j’ 


( h 0 ) 

0 /2 j’ 


and L\y — 


J 2 0 ) 
0 -/2 ]■ 


Proof Normal forms (i) are given in |5l. For (ii). Let W be a 2-dimensional symplectic subspace 
on which L has the pair of eigenvalues ±i and S(W) = W. It is known in Hamiltonian context 
that L and / can take the same normal form on W taking into account multiplication of time by 
a scalar. Equivariance property yields SL = LS. On W, L and / take the same form which gives 
SJ = JS which contradicts the fact that S is acting anti-symplectically. Thus, the minimal invariant 
subspace is four dimensional and is given by 1^ = W®W', W' = S(W). The anti-symplectic property 
implies J\w’ = -J\w while equivariance gives L\w’ = L\w = J\w- Therefore, normal forms given in 
\m show 


S|y = 


' 0 

, h 


h ^ 
0 


J\v = 


J 2 

0 


0 

-J 2 


and L\y = 


J 2 

0 


0 

J 2 


To get the same formulas for / and L given in (i) apply the change of coordinates on given by 


Zi = Wi,Z2 = W2- 

In these new coordinates S,/ and L takes the forms given in (ii). 


□ 


Note that with these conventions, L and J take the same form in both cases, and the quadratic 
part H 2 of the Hamiltonian in both is given by 

H2{Zi, Z2) ^\Zi\^ -\Z2?-, 

that is, H has a 1:-1 resonance. The higher order terms will differ for the two cases, as we see below. 
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4. Liapunov-Schmidt reduction 

The classical approach to finding periodic orbits in Hamiltonian systems is to solve a varia¬ 
tional equation on the loop space. This equation is of infinite dimension and can be reduced by 
Liapunov-Schmidt Reduction, in this section we will give an overview of that method and how to 
use it in finding periodic orbits near an equilibrium point in a reversible equivariant Hamilton¬ 
ian system. We chose the reversible equivariant case to cover all symmetry cases discussed in this 
paper. We will follow the settings given in |5) and (7) . 

Consider the vector field / : D?” — R”, which has an equivariant reversing symmetry group G. 
This implies the existence of a representation p ; G — 0{n) and a reversing sign a : G — {±1} such 
that 

fp{g) = a{g)p{g)f,ygeG. 

in the following we give briefly the main steps of the Liapunov-Schmidt reduction and details 
can be found in O . 

4.1. Defining the operator <5. Let O; be given by 

<l)(a,T)-(l + T)4^-/(u) (4.1) 

at 

where is the Banach space of IR”-valued continuous 27r-periodic functions and is the space 
of '^ 2 ji functions that are continuously differentiable, it is readily seen that zeros of O are periodic 
solutions of the dynamical system generated by / with period Now we can define the group 
action on the loop space as follows 

T:Gx^2n^^2n 

{TgU)it) = p{Y)iu{aij)t + d)), 

where g = ( 7 , 6 ) is an element of G = G k S^. Straightforward calculations imply that the operator 
O is G- reversible equivariant, that is 

<S>{TgU,T) = a{j)Tg<S>{u,T), Vg = {j,d) e G. 

The linear part of <5 is defined by 

££' - {dO)(o,o). 

It is readily verified that ^ is G-reversible equivariant. 

4.2. The splittings. Consider the splittings 

’^ 2 n - ® (ker if')-'" and‘i^ 2 ;i: - (range ® ranged?, (4.2) 

where the complements are taken with respect to the inner product 

[u,v]= (TgU,TgV)dp, 

J G 

where p is a normalized Haar measure for G and {u,v) = {u{t)Yv{t)dt. The splittings 14.21 are 
Tg-invariant. Now we define the projections 

E : ‘^^271 ^ ranged? 

I-E: ^ 2 n (rangeif')-^. 

Invariance of 14.21 under Tg implies that the projections E and I - E commute with Tg. We start 
this step by solving the equation 


E^iv+ w,t) = 0, 


6 


R. ALOMAIR AND J. MONTALDI 


for w by the implicit function theorem, where u-v+iu,vE ker ^, w e (ker The solution Vk = 
W[v,t) commutes with Tg. Thus, the Liapunov-Schmidt method reduces the original problem to 
the problem of finding the zeros of the bifurcation map which is defined by 

cp: kev£^ x |R — (rangeif)-*- 
(p{u,T] = [I - E)<^[v + W[v,j],t]. 

An important property of the bifurcation map ^ is G reversing-equivariance property, i.e 

(piTgU,j) = aiY)Tg(piu,T),yg e G. 

The last feature to be considered is the Hamiltonian structure of the bifurcation map. Using the 
implicit Hamiltonian constrain given in (7| and |5| one can show that $ is a parameter dependent 
Hamiltonian vector field. 

According to the actions of G being (anti-)symplectic we define the symplectic sign x by the 
homomorphism j : G — {±1} such that 

a){YX,Yy) = x(Y')(^(^’ylr^G. 

Therefore, the weak symplectic form Q will satisfy 

CL{gu,gv) = x{ 7 )Q(u, f),g = (r,0) e G, 

and the Hamiltonian sign is given by 

je{gu,gv) = o{.Y)xW’^^u,v). (4.3) 

In all cases we discuss kerif is finite dimensional and thus kerS£ - kerS£* and so by (71 Theo¬ 
rem 6.2] the bifurcation equation is a Hamiltonian vector field. Its corresponding Hamiltonian h 
satisfies the (semi-)invariance properties given in 14.31 restricted to kery£ i.e. 

fi(gM) = cr( 7 )X( 7 )fi(w), u e kerif, (4.4) 

where as before g = ( 7 ,0) for some 0 e S^. In practice, the function h can be computed to any finite 
degree by using normal form transformations, as described for example in (TH) (the discussion 
there is for symplectic symmetries, but is equally valid for all four cases listed in Table 1). 

5. Symplectic time-reversing involution 

In this section we prove the existence of symmetric and non-symmetric periodic solutions in a 
Hamiltonian systems with a reversing involutory symmetry acting symplectically (type SR in Table 
1). The problem was first studied by C. Buzzi and J. Lamb (5), but there is a minor sign error in 
the calculations in Lemma 6.4 which effects the statement in their Theorem 6.1. They (correctly) 
prove the existence of a three dimensional conical subspace of symmetric periodic solutions in 
a neighbourhood of the origin. Also, they find that the origin is contained in two 2-dimensional 
manifolds each containing a non-symmetric family of periodic solutions with period close to 2n. 
Using our expressions for the (semi-)invariants, we first recover their result on the symmetric solu¬ 
tions, and then we correct their Theorem 6.1 to show that generically there may or may not be two 
families of non-symmetric periodic orbits in a neighbourhood of the equilibrium point 0 depend¬ 
ing on the coefficients of the Hamiltonian. Buzzi and Lamb also distinguish between two cases, 
called elliptic and hyperbolic, distinguishing between the possibilities of the period function on 
the 3-dimensional family being monotonic or not. It turns out that this distinction coincides with 
the two cases of existence or non-existence of non-symmetric periodic orbits. 
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By the normal forms given in Lemma ini fil. we have dimker^ = 4, so we can write keri? = C^. 
Therefore, the bifurcation map is given hy 

^X IR 

(p^ZJVzh 


with Hamiltonian function 

h : X IR ^ IR, 

which satisfies Il4.4h Denote by the cyclic group generated by R, which together with gives 
XI Z 2 . The reversing symmetry R acts on by 

RiZi,Z 2 ) = {Z 2 ,Zi) 

while the action is defined by 

d{zi,Z2) = (e'®zi,e“'®Z2)- 

Let S be the ring of invariants, then one can write 

S ® S— , 

where consists of invariants and S- consists of anti-invariants. 

Lemma 5.1. LetS^ xi Z^ act on as above, then 

(1) S is the ring generated by A,B,C,D where A = |zi|^,B = |z 2 l^> C-t iD = 2ziZ2. 

(2) is the subring ofS generated by N,C,D whereN= \zif + \z 2 f,andS- is the module over 
&+ generated by the function 5 = |zip - |z 2 p- 

(3) The orbit map O: —>• IR^ defined by [zi , Z 2 ) —^ C, D) has image 

{{N, C,D]\N^>C^ + D^}. 

Note that the functions N,C,D and 6 satisfy the identity 5^ = -C^ -D^. 


The proof of this lemma is by standard algebraic computations, similar to those found for ex¬ 
ample in (9) . 

Now we can apply Lemma ET] to our Hamiltonian. The function h is -invariant, R anti¬ 
invariant and real valued. This implies there is a smooth function g satisfying 

h{zi,Z2,r] = 5 g{N,C,D,T]. (5.1) 


In order to find the periodic solutions we need to solve the bifurcation equation first. The bifur¬ 
cation equation is given by 


This can be written as 



M zh = 0 

dh 

0 dg 

dzi 


dh 

0 dg 

dZ2 



(5.2) 


We now consider, in turn, the symmetric and non-symmetric periodic orbits. 
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5.1. Symmetric Periodic Orbits. In finding symmetric periodic orbits we recover the result in (5| . 

Theorem 5.2 (Buzzi & Lamb |5]). Consider a symmetric equilibrium 0 of a reversible Hamiltonian 
vector field f with the reversing involution acting symplectically. Suppose thatDfiO) has two purely 
imaginary pairs of eigenvalues ±i with no other eigenvalues of the form ±ki, keZ. Then, the equi¬ 
librium is contained in a three-dimensional flow invariant conical subspace, given by the equation 
5 = 0, and generically this consists of a two-parameter family of symmetric periodic solutions whose 
period tends to 2n as they approach the equilibrium. 


Proof Since the Hamiltonian is R anti-invariant then all symmetric solutions are zeros of the bi¬ 
furcation equations that lie in the level set h = 0. For symmetric solutions we have 5 = 0. Therefore 
the bifurcation equation calculated in Fixi? = {(z, z) | z e C} will take the form 


Zg(Z,T) = 0 


Non-zero solutions yield g(z, t) = 0. By the formula of the reduced Hamiltonian II5.1L the lowest 
order term of the variable t is given by 

h = (|zi|^-|z 2 l^)^-t h.o.r. 


dg 1 

This implies that — (0,0) = - 0. By the implicit function theorem for each small non-zero z there 


exists a t such that (z, z) lies in a periodic orbit with period By reversing property each R sym¬ 

metric solution intersects FrxJ? in two points. Since the conical subspace 5 = 0 is 3 dimensional 
and all points in FbcR are solutions of the bifurcation equation we conclude that the conical sub¬ 


space completely consists of these periodic solutions with period close to 2n as they approach the 
origin. □ 


5.2. Non-Symmetric Periodic Orbits. We prove the existence of two families of non-symmetric 
periodic solutions under suitable conditions on the coefficients of the Hamiltonian. This result is 
fairly different to the one in |5|. To prove the existence of non-symmetric solutions one needs to 
solve the bifurcation equation without any symmetry conditions. By calculating the partial deriva¬ 
tives of g the bifurcation equation will be 


— = Zi{g-t5giv) + Z25(gC-2gD) ^0 
OZi 

dh 

■^^Z2{-g + SgN) + Zi5{gc-lgD)^0 
dZ2 


dg dg dg 

where Sn = gc = -^ and go = ^- Multiplying the first equation by Zi 
by Z 2 we get 


and the second one 


l■Zll^(g + <5giv) + ■^lZ25(gc- ign) = 0 (5.3) 

Iz 2 l^(-g + 5giv) + ZiZ25(gc - igo) = 0 (5.4) 

By adding II5.31 and 115.41 we have 


5(g -t NgN -t (C -t iD]{gc - igo)) ^ 0 


(5.5) 


Taking the imaginary part of the above equation gives 


^gc - CgD = 0 
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and when C,D^ Owe can write that equation as 


gc _ gD 

C ~ D 

Therefore equation <5.5< will be 

(5.6) 

5{g + NgN + Cgc + Dgo) = 0 

By subtracting <5.4< from <5.311 we have 

(5.7) 

Ng + 6^gN^0 

this can also be written by the formula 

g _ gN 

(5.8) 

52 “ N 

Substituting <5.6< and <5.9ll in <5.7< yields 

gN _ gc 

(5.9) 

N ~ C 

Thus 

g _ gN _ gc _ gD 

(5.10) 

52 “ N ~ C ^ D 

which is equivalent to 

N C D 

(5.11) 

gN gc gD ’ 

(5.12) 

In order to prove the existence of non-symmetric periodic solutions to the original Hamiltonian 
system we need to prove the following lemma. Let 

g]v( 0 ) = n,gc( 0 ) = c,gi)( 0 ) = d. 

Lemma 5.3. Ifn,c and d are not all zero then there exists a unique solution in IR^ = 
for the system of equations 

{t,N, C, D) -space 

g + NgN + Cgc + Dgo = 0 

(5.13) 

Ngc + CgN = 0 

(5.14) 

Dgc - CgD = 0 

(5.15) 

NgD + DgN = 0. 

(5.16) 


Proof. It is clear that the last three equations are not independent hut we will use them all to make 
up for the special cases when one of the numbers n,c or d is equal to zero. Suppose that n ^ 0. 
Then we only need to solve <5.131 1 , <5.141 1 and <5.1611 . In order to apply the implicit function theorem 
we need to study the following Jacohian matrix with respect to t, C, D and N 


f 1 

2 

c 

d 

n A 

0 

n 

0 

c U( X| 

0 

0 

n 

d ) 


Since n ^ 0 then the matrix X is non-singular. Therefore hy the implicit function theorem there 
exists a unique curve S = SiN), with dS[0) = -X~^Y, that solves the system. If n = 0 hut c 5 ^ 0 
we can choose equations <5.13< . <5.14< and <5.15< . Solving hy the implicit function theorem gives a 
unique solution S = S(C). A similar argument can he used for the remaining cases. 

□ 


Now we state and prove the main theorem about the existence of non-symmetric periodic solu¬ 
tions for the given reversible Hamiltonian system. 
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Theorem 5.4. Suppose that ^ + (fi, then there exist the symmetric Liapunov centre families of 

periodic solutions filling the set 5 = 0 described before. Moreover, 

i) Ifn^>c^ + d^ then there exists two families of non-symmetric periodic orbits for the Hamil¬ 
tonian system distinguished by the sign of 6. The period of the periodic solutions converges 
to 2n as the solutions tend to the origin. 

ii) Ifn^<c^ + d^ then the only periodic orbits with period close to 2n in a neighbourhood of the 
origin are the symmetric ones. 

Proof. To prove the existence of non-symmetric periodic orbits we have to solve the equations 
II5.13MI5.14I I. II5.15I I and 15.1611 . By the condition n^ ^ c^ + d? we have that n, c and d cannot all be 
zero. Applying Lemma lS^l we have a unique solution for those equations. Therefore, we can write 


N 

gN 


— --—-t 
gc go 


(5.17) 


which is equivalent to N = gj\ft,C = -get and D = -got. To get non-symmetric solutions we 
should have 5^ = -C^ - > 0. This implies 

igN^ - gc^ - gD^) > 0, for t 0, 

and therefore, - gc^ - gn^ > b. Taking the limit at the origin gives n^'^ c^+d^. We conclude that 
non-symmetric solutions exist when n^ > c^ -\- d^ and split into two families according to 5 being 
positive or negative. On the other hand, when n^ < c^ d^ the only periodic orbits with period 
close to 2n in a neighbourhood of the origin are the symmetric ones. □ 

5.3. Period Distribution within the Famiiy of Symmetric Periodic Solutions. Following the ar¬ 
gument given in Buzz! and Lamb O , we describe the structure of period distribution for symmet¬ 
ric periodic solutions. According to FbcR being two dimensional the level sets of the period will 
be given by t = t(x, y). If we change the coordinates in a neighbourhood of the origin such that 
T = -I- £ 2 y^ with Cj - ±l,where the sign depends on the details of h and H. One can give the 

following definition: 

Definition 5.5. The level sets of the period t can be of two types: 

(1) elliptic vdoen £i £2 = 1. In that case the level sets of the period form circles and t increases 
or decreases monotonicaily with increasing radius. 

(2) hyperbolic when £i £2 = -1. Here the level sets of the period form two families of hyperbo¬ 
lae, one family with positive increasing t and one with negative decreasing t. 

Now we can prove the following proposition: 

Proposition 5.6. Dependingon the quartic terms of the Hamiltonian function (and quadratic terms 
of the function g), among the three dimensional surface of symmetric periodic solutions near the 
equilibrium point, then the level sets ofr are elliptic when n^ > c^ -t d^ or hyperbolic when n^ < 
c^ + d^. 

Proof. As discussed in the proof of the existence of symmetric periodic solutions, t(x, y) can be 
calculated using the equation g{z,T] = 0, with z = x + iy. Using our variables N,C and D and 
depending on the quadratic terms of that equation we have 

g[N,C,D,T) = 0 

T 

nN cC -t dD _= — 

2 

2n(x^ + y^) -)-2c(x^ -y^) -4d(xy) = 
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By the Morse Lemma the shape of t[x, y) near the origin is given hy the determinant 

D = 4^{n^-c^-d^). 

Therefore the family of periodic orhits is elliptic when {n^ - - d^) > 0 or hyperbolic when - 

- d^) < 0 . □ 

Accordingly, one can easily deduce the following corollary. 

Corollary 5.7. The two dimensional families of non-symmetric periodic orbits given in Theorem \5^ 
exists if and only if the three dimensional family of symmetric periodic orbits is of elliptic type. 

6. Anti-symplectic involution 

In this section we analyse the problem of existence of periodic orbits in a Hamiltonian system 
which is equivariant under the action of an anti-symplectic involution S (type AE in Table 1). This 
was studied by J. Li and Y. Shi in (TS), but that paper contains a number of errors. Firstly, the 
form of the Hamiltonian is not sufficiently general, for example the polynomial function h = DN 
satisfies the symmetry of the problem but is not in the form assumed in (TSl . This effects the results 
significantly and the general form of the Hamiltonian makes the calculations more difficult. There 
is also a serious error in the proof of their Lemma 5.3. As a result we consider the problem anew. We 
use a different basis from (13] , so the invariants and anti-invariants are different, and we determine 
a general formula for the reduced Hamiltonian. Firstly, we find that no symmetric periodic orbits 
can occur generically (opposite to the result claimed in (131). Secondly, we prove the existence of at 
least two and at most 12 families of non-symmetric periodic solutions near the equilibrium point. 

An immediate consequence of our assumptions is the Hamiltonian being S anti-invariant (as 
pointed out in Table 1). By the normal forms given in Lemma IXTl ii) we have dimkerif = 4 i.e. 
kerif = C^. The bifurcation equation is given by the formula 

^X IR ^ 

(p = 2jyzh 

with the Hamiltonian 

h:C^ xU^U, 

where / is the structure map. Now let’s define the actions of Z| x on by 

S(Zi,Z2) ^ iZ2,Zi] 
d{zi,Z2) = {e‘^zi,e~‘^Z2) 

Now we study the set of (anti-)invariants and find the appropriate formula for h. 

Lemma 6.1. ForS^ x acting on as above, then 

(1) TheS^ X Z 2 invariant functions are generated by N,C,D^ where 

N = \zi\^ + \Z2\^ ,C + iD = 2 ziZ2. 

(2) theS^ invariant but Z 2 anti-invariant functions are generated by 6 , D where 

5 ^ lZi\^-lZ2f. 

According to that the Hamiltonian h will take the form 

h = dg^ [N, C,d\t) + Dg^[N, C,D^,t). 
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The bifurcation equation will be given by 


dh 

dzi 


= zig^ +5^-iz2g^ + D^ = Q 
OZi OZi 


dh 

dZ 2 


dZ2 


2 ^6 

iZig^ + D— ^0 
dZ2 


( 6 . 1 ) 

( 6 . 2 ) 


6.1. Symmetric Periodic Orbits. Symmetric periodic solutions of that equivariant Hamiltonian 
system lie in the set FixS = {(z,z),z e C}. Moreover, by anti-invariance, that is ho S = -h, all sym¬ 
metric solutions will be in the level set b = 0 . In order to get the symmetric periodic solutions we 
need to solve the bifurcation equation calculated in FixS. Consequently, one needs to solve I l 6 .lt 
and 16.21 with conditions: 6 = D = 0 and N= C. Thus, 

zig^ - iz 2 g^ = 0 
-Z 2 g^-izig^ = Q 

By multiplying 16.311 by Zi and 16.411 by Z 2 we get 

|Zi|^g^-iZiZ 2 g^ = 

-|Z 2 l^g^-iZiZ 2 g^ = 

Adding and subtracting these two equations yields 

5g^-i(C+iD)g2-0 

Ng^=0. 

With the conditions 5 = D = 0 we have 

Cg ^=0 

Ng^=0. 

Since we are looking for nonzero solutions then N = C 7 ^ 0 and therefore, solutions are common 
zeros of g^ and g^ in a neighbourhood of the origin. But g^ and g^ are independent functions and 
generically the only common zero in a neighbourhood of the origin is 0 itself. As a result there are 
no symmetric periodic orbits for the given Hamiltonian system. 

Remark 6.1. Another way to see the non-existence of symmetric solutions in that system is by 
using a Liapunov function. Consider the Hamiltonian given by the formula H = 6 [ai + biN+ciC + 
■ ■ ■) + D{a 2 + b 2 N + C 2 C -\—). Restricting the Hamiltonian system on the two dimensional invariant 
space FixS gives 

X = 2y (rJi -H 2.[bi -t Ci)[x^ -t y^) + • • •) -t 2x (£12 "t 2(^2 + Cz)T^) + ■ ■ ■) 

y = — 2 x (rJi + 2 (hi -1- Cj) (x^ -1- y^) -t • • • ] -1- 2 y (£J2 “t 2(^2 + C2) (x^ + y^) + —) 

Easy computations show that the eigenvalues of the linear system are A = 2 (fl 2 ± i)- In order to 

get periodic orbits we should have fl 2 - 0 and the system would be written as 

X — 2 y [rjj + 2 ,(b\ + ci) (x^ + y^) -1- ■ ■ ■) -t- 2 x [2(^2 + C2)(x^ + y^) + —) 
y — — 2 x [ni + 2 ,[b\ + ci)(x^ + y^) -1- ■ ■ ■ ] -t- 2 y [2(^2 + C2)(x^ + y^) -t- ■ ■ ■ ]. 



(6.3) 


(6.4) 

0 

(6.5) 

0 

(6.6) 
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Consider as Liapunov function V = + y^. Differentiating V in the direction of the Hamiltonian 

vector field yields 

V = 2xx + 2yy 
= 8[x^ + y^fib2 + C2). 

The number ha + Ca is generically non-zero and therefore V is non-zero. This means the sign of V 
(either positive or negative) is constant along any trajectory, so that the trajectory cannot be closed. 
Thus, the system does not have any symmetric periodic orbits. 

6.2. Non-symmetric Periodic Orbits. For this case we only need to solve the pair 16.11 1 and 16.21 
without any extra conditions. Multiplying 16.111 by Zi and 16.21 1 by Za gives 

|Z 1 -t 5 -t gcZiZ2 + gjj22D[-iziZ2)) - iziZ2g^ 

+ D[gff\Zif + glziZ2 + g^22D{-iziZ2)] ^0 (6.7) 


- \z2fg^+S[glj\z2\^ + gcZiZ2 + g^22D(-iziZa)) - iziZ2g^+ 

+ D[gff\z2f + glziZ2 + g^22D{-iziZ2)] ^0 ( 6 . 8 ) 

By adding theses two equations we have 
d [g^ + Nglj + (C + iD)g^ + 2g^2 {-iD)[C + iD)] - i (C + iD)g^ 


+ D[Ngff + (C-E iD]gQ + 2g^2i^iD)iC + iD)] — 0 (6.9) 

The real and imaginary parts of equation 16.91 are 

5 (gi + ATg^ -E Cgi + 2g^2 o2) + Dg2 + D (iVg^ -E Cgl + 2g2, D^) = 0 ( 6 .10) 

5 [Dg^ - 2g^, CD) - Cg2 + D [Dg2 - 2g2, CD) - 0 (6.11) 

The last equation to be considered comes from subtracting 16.81 from 16.71 and it will take the form 

A^gi+<5"g^ + Ddg ^-0 ( 6 . 12 ) 

This means hnding non-symmetric solutions of the Hamiltonian system will be by solving the 
triple 16.101 . 16.111 and 16.121 . Clearly the system is singular at the origin and can be studied us¬ 
ing a blow-up method. For that purpose define the new coordinates [u, v, w, t, x) by 

N=rv, C=ru, D-rw, 


T = rt, 6 = rx, 

combined together by the relation + x^ according to the relation = 6^ + C^ + D^. 

Substituting these new coordinates in 16.101 . 16.111 and l6.121 gives 

r (i;g^ -E rx^glj + rxwgjj] = 0 

r (x(g^ + r vglf +rugc + 2r^ w^gD 2 ] + w[g^ + r vg% +rugl + 2r^ w^goz)] = 0 

r [xig^rw - 2 r^ujugjj 2 ) - ug^ + w[rwg^- 2 r^uwg^ 2 )] = 0 


(6.13) 
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We are interested in the non-zero solutions, i.e. r ^ 0. The first step is to divide by the common 
power of r in these equations and the second step is to apply the implicit function theorem. For 
simplicity we can write the Taylor series for the functions and as 

g^ = — -t ill Ai CiCdi -[-••• 
g^ = ^2^ + 0-2^ C 2 C + d2D^ -!-••• 


which with the new coordinates take the form 

- rg^ = r(^-i-air’-tCiU-t(iirm^H—) 
g^ = rg^ = rib2t + a2V + C 2 U + dzr + ■ ■ ■) 


Accordingly, the system US. 131 will be written as 

r^[vg^+ x^gl + xw£] = 0 

(x(g^ -h vgl + ugl + 2 rw^gl^ 2 ) + + vgl + ugl + 2r - 0 

[xiwgl - 2 rwugl^ 2 ) - uf + wiwgl - 2 ruwgl^ 2 )) = 0 
Note here that g^ = gj etc. Dividing by and substituting r = 0 yields 

t 2 

v{- + aiv + ciu) + aix +a2xw = 0 


(6.14) 


x{- +2aiv + 2c\u) + w[h 2 t + 2 a 2 V + 2 c 2 u) = 0 

C\xw- u{b 2 t + a 2 V + C 2 U) + C 2 W^ = 0. (6.15) 

Clearly, the system can not be solved by the implicit function theorem at this point in the argu¬ 
ment. As a result we will use a different technique as illustrated in the next section. We will show 
that 16.151 has non-degenerate solutions, then apply a continuation argument to show 16.141 has 
solutions when r > 0. Adding the relation between the variables N, C, D and 5 gives us the system 


t 2 

vi-+ aiv + Ciu] + aix +a2XW-0 

t 

x[- + 2 aiv + 2 ciu) + w{b2t + 2a2V + 2 c2u) = 0 
CiXU> - U[b2t + a2V + C 2 U) + C2lU^ - 0 


2 2 2 
u + w; + X ■ 


v^ = 0. 


(6.16) 

First of all we want to count the number of all solutions of the system 16.161 . For that purpose 
we need the following theorem. 


Theorem 6.2 (Bezout’s theorem). Suppose n homogeneous polynomials onC in n + \ variables, 
of degrees di,d 2 ,..,dn, that define n hypersurfaces in the projective space of dimension n. If the 
number of intersection points of the hypersurfaces is finite, then this number is did 2 -dn if the points 
are counted with their multiplicity. 

For more details and proof see for example (10) . 

The system 16.161 consists of four homogeneous equations each of degree two with five vari¬ 
ables. According to Bezout’s Theorem we have 16 complex solutions for that system and can divide 
them into two main types: solutions when v = 0 and solutions when v 
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6.2.1. Solutions when v = 0. In that case, algebric calculations give a total of three different solu¬ 
tions: 


( 1 ) {t e K, n = 0 , ui = 0 ,x = 0 }, 

(2) {t =+2i^w, u = ±iw, w eU, X = 0}. 

Now we want to study the multiplicity of each solution. Consider the Jacohian matrix for the sys¬ 
tem II6.16> with respect to v, t, u, w, x 

' jt+2aiv+ciu CiV azx Zaix+Uzw 

J— 2aiX+2a2W ^x+b2W 2ciX+2c2W 2a2V+b2t+2c2U ^t+2aiv+2ciu 
-a2U -b2U -a2V-b2t-2c2U Cix+2C2W Ciw 

< -2v 0 2u 2w 2x 


Substituting the values of the first solution and the condition h = 0 in the Jacobian matrix yields 


2 

0 

0 

0 

0 J 

0 

0 

0 

b2t 


0 

0 

-b2t 

0 

0 

lo 

0 

0 

0 

0 J 


To get the appropriate square submatrix we eliminate the second column because t is non-zero 
and get 

(^t 0 0 0 \ 

0 0 b2t \t 

0 -h2t 0 0 

,0 0 0 0 ; 

This matrix is of rank three and therefore this first solution is not simple. To study its multiplicity we 
need to study the behaviour of system I l6.16t near a solution point for example say [v, t, u, w, x) = 
(0,2,0,0,0). Consider the system 

ll(l + ayV + CiU] + ayX^ -1- a2XW - Cl 


x(l -i- 2 ai v + 2CiU] + w[ 2 b 2 + 2 a 2 V + 2 c 2 U) = £2 

Cixw- u{2b2 + a2V + C 2 U) + C2iv^ = £3 

+ x^ - (6.17) 

Near the point [v, t, u, w,x) = (0,2,0,0,0) the first equation can be solved by the implicit function 
theorem for v, the second for x and the third equation for u. As a result we end up with solving the 
equation 

+ f{w) = £4, 

where f{w) is a function constructed by substituting the solutions from the implicit function the¬ 
orem in equation <6.1711 . Clearly / is of order greater than one. So, the least order coefficient is 
and so the studied solution is of multiplicity two. 


Regarding the multiplicity of the second and third solution we should assume that w ^0 for a 
non-zero solution; for simplicity let le = 1. The Jacobian matrix will take the form 


'+C2i/b2 ± Cii 

0 

0 

0 

02 

2fl2 

h2 

2C2 

0 

+C2ilb2±2cii 

+a2i 

+b2i 

0 

2C2 

Cl 

0 

0 

+2/ 

2 

0 j 


y 
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Since tc = 1, we can omit the w- column and get 


h - 


'+C 2 i/b 2 ± Cli 

0 

0 

2 a 2 

b 2 

2 C 2 

+ CI 2 i 

+b 2 i 

0 

1 0 

0 

+ 2 i 


a2 ^ 
+C2ilb2±2cii 
C\ 

0 


det /2 = + b\c\- 2 b 2 CiC 2 + cf)/h 2 . 


We can assume that this result is non-zero and therefore the second and the third solutions are 
simple. We conclude that the case v = 0 corresponds to four solutions where the first solution is 
doubled hut the others are of multiplicity one. Note that v = 0 implies N = |zi|^ + \z 2 f = 0. Thus, 
these four solutions won’t he counted as periodic solutions of the given system, hut will help us 
find out how many non-zero periodic solutions there are. 

6.2.2. Solutions when v^O. There remain 12 solutions for the case v 7 ^ 0 according to Bezout’s 
theorem . The following proposition guarantees a minimum of two real solutions for the system 
EEl. 

Proposition 6.3. For any choice of coefficients {ai, 02 , b 2 , ci, C 2 } the two points 

{n G K*, f - -Aa\V, u= w -0,x = ±v}, 
satisfy equation 16.161 , when vt^o. 

Proof Straightforward calculations yield the result. □ 

In order to find out more about the maximum number of real solutions we can find we will use 
a numerical approach. We choose various values for the constants in the system 16.161 and then 
solve the equations using Maple. Since we are interested in solutions with y 7 ^ 0, we put y = 1 for 
simplicity. These numerical calculations suggest that the system can have a maximum of eight 
real solutions, including the two analytic solutions given by Proposition l6.3l In addition, there are 
examples of systems with four or six real solutions. Our aim is to prove that for each of these cases, 
the solutions are non-degenerate. Then, under any perturbation of the set of coefficients there still 
exist (nearby) real solutions (i.e. periodic solutions). In the following we study an example of each 
set of coefficients that has two, four, six or eight real solutions for the studied system 16.161 . Then, 
we check their non-degeneracy conditions. Note that ail numbers are rounded to four decimal 
digits. 

Example 6.4 (A system with two real solutions). Consider the set 

R- {ai- 1, fl 2 -5,b2 = 1, Cl = 2, C 2 = 2, f = 1}. 

The corresponding system has only two real solutions 

{t- -4, u-0,w = 0,x-±v-±l}, 

which are those given in Proposition 16.31 The remaining 10 solutions are non-real. In order to 
check the non-degeneracy condition, we need to study the proper submatrix of / for each solu¬ 
tion and ensure that its determinant is non-zero. Substituting the values given in R and the two 
solutions in / yields 


/i- 


0 

0.5 

2 

±5 

+21 

+2 

±0.5 

±4 

6 

0 

0 

0 

-1 

±2 

0 

1-2 

0 

0 

0 

± 2 ; 
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(° 

2 

±5 

±2^ 

±2 

±4 

6 

0 

0 

-1 

±2 

0 

1-2 

0 

0 

±2; 

det/11 = 

±20. 



Since i 7^ 0, we omit the t-column and we have the suhmatrix 


/ii = 


Therefore, these two solutions are non-degenerate. 

A similar argument is used in the remaining examples to prove the non-degeneracy of solutions 
in each case. 

Example 6.5 (A system with four real solutions). Let the set of coefficients in the system <6.161 1 he 

R = {a-i = 1, a 2 = 5,b2- -2, ci = 2, C2 = 2, h = 1}. 

the associated system has 8 non-real solutions and only four real solutions and the real ones are : 

(1) {t = -A:,u = 0,w = 0,x = ±v = ±1] 

(2) {t = 1.5602, u = -0.9681, w = +0.0855, x ^ +0.2354} 

Suhstituting R and the first two solutions in the matrix / we get 


/i- 


Now we can choose the suhmatrix /n hy omitting the second column because t is non-zero and 
we find its determinant to be det Ji 1 = +692 0. In the same way we can study the third and fourth 

solutions to get 

{ 0.8439 2 +1.1772 ±0.8985^ 

±1.3262 +1.2839 3.0071 -1.0924 

4.8406 1.9929 +0.8130 ±0.1711 

V -2 -1.9362 ±0.1711 ±0.4709/ 

We have det/22 = ±35.6351 5^ 0. 

Since the determinants are non-zero, all four solutions are non-degenerate and we can find an 
open set of coefficients that give four real solutions. 


( ^ 

0.5 

2 

±5 

±21 

±2 

±0.5 

±4 

18 

0 

0 

0 

-13 

±2 

0 

1-2 

0 

0 

0 

±2, 


/22 - 


Example 6.6 (A system with six real solutions). Let 

R = {ai- -2, fl2 - -11.^2 = -5, Cl = 1,C2 = 2, c = 1}. 

The system <6.161 with those coefficients has only six real solutions: 

(1) {t = ^,u = Q,iv = 0,x = ±v = ±1] 

(2) {t = -2.5592, u = 0.0346, w = ±0.4980, x ^ +0.8665} 

(3) {r ^ -3.7663, u = 0.1529, w = ±0.8984, x ^ +0.4118} 

The non-degeneracy of the above solutions can be studied in pairs. Firstly, we study the determi¬ 
nant of the appropriate matrix Ju associated to the first and second solutions. 


/ii = 


0 

1 

+ 11 

+41 

+4 

±2 

-62 

0 

0 

51 

±1 

0 

1-2 

0 

0 

±2; 
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det/ii = +20816. 

Similarly, for the rest of solutions we have 


1-5.2450 1 +9.5314 

+7.4900 +0.2590 -9.0658 

0.3804 -1.9342 +1.1255 

, -2 0.0692 +0.9960 


+ 2.01201 
-5.2104 
±0.4980 
+ 1.7330, 


det/ 22 ^+164.8123 


/33 


-5.7303 1 +4.5299 

+ 18.1165 ±2.7698 -2.5567 

1.6818 -8.4433 ±3.1816 

, -2 0.3058 ±1.7967 


+8.23461 

-5.5774 

±0.8984 

+0.8236, 


det /33 ^ ±1827.2294. 
As a result all real solutions of this case are non-degenerate 


We end with an example of a system with eight real solutions, which is the largest number of 
real solutions we found using numerical calculations. 


Example 6.7 (A system with eight real solutions). Let 

R-{ai- l,a2 = -4,^2 = -1, Ci = 1, C 2 = 2, f = 1} 
The corresponding real solutions are only eight and they are 

(1) {t =-4, a = 0, lu = 0,x = ± 1 '= ±1} 

(2) -4.9432, u = -0.2615, w = ±0.2274, x ^ +0.9380} 

(3) -2.8537, u = 0.8527, w = ±0.4155, x ^ ±0.3165} 

(4) -6.4260, u = 0.2940, w = ±0.8063, x ^ +0.5133}. 

The non-degeneracy conditions are 



( ® ^ 

+4 ±21 


/ll = 

±2 ±2 

0 

1 


0 0 

±1 0 



0 

CM 

1 

0 ±2J 



det /11 ^ 

±4 


7-0.7331 

1 

±3.7521 

+2.78571 

+3.6953 

+0.9664 

-4.1029 

-0.9947 

-1.0461 

0.1029 

+0.0284 

±0.2274 

-2 

-0.5231 

±0.4548 

+ 1.8760, 

det J22 — +13.8083 


1.4259 

1 

+ 1.2659 

+ 1.0293^ 

+2.6915 

±2.2951 

-1.7353 

2.2786 

3.4110 

-2.2647 

±1.9787 

±0.4155 

-2 

1.7055 

±0.8311 

±0.6329, 


det J33 ^ +43.7450 
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^-0.9190 1 

+7.4767 +2.1984 

1.1761 -3.6021 

, -2 0.5881 


+2.0533 +4.2517^ 
-0.3979 -0.6249 
+2.7117 ±0.8063 
+ 1.6125 +1.0266, 


det/ 44 ^ ±111.6657 

Therefore, all eight solutions are non-degenerate. 


6.3. Conclusion. Bezout’s theorem guaranteed a total of 12 solutions for the case v ^ 0, hut nu¬ 
merical calculations found at most eight of them to he real (and at least two). The last thing to 
consider is the effect of the addition of higher order terms to the system 16.161 when solving hy the 
implicit function theorem. We will choose one of the previous examples and prove the existence 
of periodic orhits in that system and the rest can be done in the same way. 

We consider the solution point [t, u, w, v,x, r) = (-4,0,0,1,1,0) as a candidate. We want to ap¬ 
ply the implicit function theorem on the system in a neighbourhood of that point. Note that the 
functions g^,g^ are given by 


g^(A^,C,D^T) = ^ + aiN + CiC + diD^ + eiN"^ + fiNC + giNr + ■ ■ ■, (6.18) 

g^{N,C,D^,T) = b2T + a2N + C 2 C + d2D^ + ezN^ + f2NC + g2NT + ■■■ (6.19) 

In our new coordinates 16.1811 and 16.191 will take the form 

g^{N, C,D^,T) = r[- + aiV + ciu + dirw^ + eirv^ + firvu + girvt -\—], ( 6 . 20 ) 

g^{N,C,D^,T) = r[b 2 t + a 2 V + C 2 U + d 2 r + e 2 r + f 2 r vu + g 2 r vt + ■ ■ ■], (6.21) 


therefore, the matrix formula associated to the implicit function theorem calculated at the point 
(-4,0,0,1,1,0) will be 


0 

CM 

5 

2 

3ei-8gi ^ 

2 

4 

6 

0 

3ei-8gi 

0 

-1 

2 

0 

0 

-2 

0 

0 

2 

0 


The matrix X is invertible and by the implicit function theorem we can solve v, u, w, x as functions 
of r. The linear part of the Taylor series of those solutions is determined by the matrix 


7/5 

-9/10 

-4/5 

-7/5 



( 3ei 

-8gi ^ 

-2/5 

2/5 

-1/5 

2/5 



3ei 

-8gi 

-1/5 

1/5 

2/5 

1/5 




0 

7/5 

-9/10 

-4/5 

-9/10 




0 


Those solutions can be written as functions of r as follows 

H(r) = 1 - i (3ei - 8 gi)r + h.o. t. 

x(r) = 1 - i (3ei - 8 gi) r + h.o. t. 
v-w = 0 . 
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Converting back to our basic coordinates N, C, D, 6 gives 

1 2 , 

N = rv = r - -(3^1 -8gi)r + h.o.t. 

1 , 

6 = rx = r- -(3ci -8gi)r + h.o.t. 

C = D = 0 

Tbis curve of solutions gives a one parameter family of periodic orbits for tbe equivariant Hamil¬ 
tonian system. Similarly one can prove tbe existence of one parameter family of periodic solutions 
for each case studied before because of tbeir non-degeneracy conditions. Accordingly, we state tbe 
following result. 

Theorem 6.8. Consider an equilibrium point 0 of a C“ equivariant Hamiltonian vector field f, 
with the the symmetry S acting anti-symplectically and = I. Assume that the linear Hamiltonian 
vector field L has two pairs of purely imaginary eigenvalues ± i and no other eigenvalues of the form 
±ki, keZ. The reduced Hamiltonian is in the form h = 5g^{N,C,D^,T) -\-Dg^{N, C,D^,t). Then 

(1) For an open dense set of coefficients (ai, fl2. i’2. ci, C2) there exists a neighbourhood ofO with 
no symmetric periodic orbits and at least two and at most 12 non-symmetric periodic solu¬ 
tions of the equivariant Hamiltonian system. 

(2) There exist open sets of coefficients Ui [i = 1,2,2,A), such that for coefficients inUi thereare 
precisely 2i non-symmetric periodic orbits of period close to 2n as they tend to zero. 

7. The combined case x 

It is natural at tbis point to ask about periodic orbits of a system possessing both tbe symmetries 
studied above. Consider now a reversible equivariant Hamiltonian system under tbe action of tbe 
group G = Z 2 X Z|, where R and S are tbe involutions defined in Section 5 and Section 6 respec¬ 
tively. In this section we prove the existence of families of periodic solutions in a neighbourhood 
of the origin in that system. 

On the reduced Hamiltonian is a special case of the Hamiltonian in Section 5 and it takes the 
form 

hizi,Z2,T) = 5giN,C,D^,j). 

Accordingly, the bifurcation equation will be 

6[g + NgM + Cgc + 2D^gjj2] = 0, 

Ng+[N^-C^-D^]gN = 0, 

5Dlgc-2CgD2]^0. (7.1) 

Now we classify the solutions according to their symmetry type. 

7.1. Periodic orbits in the conical subspace 5 = 0. Substituting 5 = 0 in the system of equations 
17.H yields 

Ng = 0. (7.2) 

For R symmetric solutions one needs to solve 17.2t in FixR. This implies 

giZ,T)=0, 

which can be solved for t = t(z),z e Fixl? by the implicit function theorem. This means any pe¬ 
riodic orbit in the subspace 5 = 0 has symmetry R. Solving equation l7.2ll for z e FixS gives one 
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Figure 1. Fixed point spaces on the torus (0i,02)- 

periodic orbit of symmetry S and it is therefore SR symmetric. Moreover, solving equationl l7.2t for 
z e Fix(S, 7 t) gives another orbit with symmetry S. 

7.2. Periodic orbits in 5 0. It remains to study the existence of SR solutions which lie in the 
subset 5 7^ 0. Clearly FhcRS = {(zi.za) I Zi,Z 2 £ D?} which implies D - 0 and therefore the system 
Il7.1i takes the form 

g + Ngjv + Cgc ^ 0, 

Ng+[N^-C^]gN^0. (7.3) 

Eliminating g from both equations gives 

C(Ngc + CgN)^0. (7.4) 

If C = 0, then by the fact (0) - | we can solve using the implicit function theorem. Now if C 5^ 
0, and g]v(0) = n,gc(0) = c are not both zero then, the svstem <7.3> can be solved by the implicit 
function theorem. By the argument used in Theorem 5.2 we conclude that SR periodic solutions 
exist when > 0. The following theorem describes the families of periodic solutions exist in 

this system. 

Theorem 7.1. Consider a symmetric equilibrium 0 of a Z 2 x reversible equivariant Hamiltonian 
vector field f where R is a reversing involution acting symplectically and S is an involution acting 
anti-symplectically. Suppose that Df(0) has two purely imaginary pairs of eigenvalues ±i with no 
other eigenvalues of the form ±ki, keZ. Also, denote gj\i(0) = n andgc(0) = c. Then, 

(1) there exists a two-parameter family ofR symmetric periodic solutions in the conical subspace 
6 = 0 with two of them having extra symmetry S. The period of all orbits tends to 2n as they 
approach the equilibrium. 

(2) there exist two Liapunov centre families of SR symmetric periodic solutions in the open sub¬ 
sets 7^0 provided that n^ - c^ > 0 —one with 8 > 0 and one with 6 <0. These two families 
are exchanged by both involutions R and S. 

Finally we illustrate the relation between fixed point spaces of the involutions R,S,SR and (S, n) 
geometrically. Buzzi and Lamb (5| show that the intersection between the cone 6 = 0 and the unit 
sphere in is a torus T parametrized by two angles ( 61 , 62 ) and draw FixR on T. In addition to 
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that we show the intersection between FixS and the torus T is given hy the line 62 = - 6 i. Also, we 
plot Fix(S,:;r) = {( 61 , 62 ) = ( 6 i,n- 6 i)} on the torus. The last thing is to intersect FixSi? with T which 
gives a total of two points (0,0) and (n, n) (shown as large dots in the figure). 
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